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There are known problems of Lorentz-Dirac equation for moving with ac- 
celeration charged particle in classical electrodynamics. The model of ex- 
tended in one dimension particle is proposed and shown that electromagnetic 
self-interaction can lead (with appropriate choice of retarded and advanced 
O I interactions) to zero change in particle momentum. The hypothesis is for- 

^ ■ mulated: all relativistic internal forces of various nature can give zero change 

Q ■ in particle momentum 

l> : 03.50.De 

^ 0. 

Qh! Since the famous Dirac's paper on relativistic radiation reaction in clas- 

^ I sical electrodynamics, many textbooks and research articles were published 

on that theme. Among them are [1-11], where one can find the discussion of 



> 



X 



the related problems: mass renormalization and its nonuniqueness, runaway 
solutions and the use of the advanced interaction. ( These problems of radi- 
^ ', ation reaction one can find also in other classical field theories - scalar field 

theory, gravitational theory, etc.) 

First of all we give brief review of some problems (sections 1-3) and then 
in section 4 consider the model of an extended particle. 



It is well known that the x-component of electric field E produced by 
point particle with charge e = 1 moving along x-axis is ( Lienard-Wiechert 
solution; we take the units c = 1) 



P {l + Nv')N 

[l - Nv'){x - R'f ^ ■ ' 
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here R = R(t') - is the trajectory of particle motion and v' = v{t') = dR' /dt' 
- is the particle velocity, t' - is the retarded time: t' = t ~ \x — R'\, N - is the 
unit vector: N ^ {x - R')/\x - R'\ 

To consider self-interaction one must expand field E near the particle. 
One can do it in different ways. For example: 
1. Let e 



t - t' (i.e. \x - R! 

1 



<E>- 



= e), e — > 0. Then 



2v 



2i>(l + 't;2) 



+ 



e2(l--t;2) e(l-'t;2)2 ' [{I - v"^)' 
2. Let e = |x - R{t)\, e ^ 0. Then 

1 



+ 



<E>^- 



+ 



2v (vy2v 



3fl 



,2^2 



+ 



^2)3 



3. Let e — \x — R'\ 



{l-Nv') 
{1+v'N) 



,6^0. Then 



:l2) 



<E>^- 



= 



[lA) 



here dot means differentiation with respect to t, all quantities in < > are 
taken at moment t. 

Infinite at e ^ terms in (f .2-4) lead to mass regularization in equation 
of particle motion. Finite at e ^ terms in (1.2-4) are what is usually called 
"the radiation force". Finite terms in (1.2-1.4) disagree each other and only 
finite term in (1.3) gives the known expression for Lorentz-Dirac radiation 
force. 

Thus different ways of expansion of particle self-field may lead to different 
"radiation forces". From mathematical point of view it is obvious result - 
the expansion near infinity is not well defined operation. 

This is a one of various aspects of the problem of mass regularization and 
its nonuniqueness (see also [2,17]). 
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Among other problems, the one, which, to my opinion, draws not enough 
attention in the hterature - is the problem of instability of preaccelerative 
"physical" solution of relativistic Lorentz-Dirac equation. 

Let the point particle with mass m and charge e move under some external 
force F along x-axis. The relativistic Lorentz-Dirac equation reads: 



mv 



(1 - (^;/c)2)3/2 



F + 



2e^ 
3c3 



+ 



3v{vy 



(l-(v/c)2)2 c2(l - (^;/c)2)3 



(2.1) 



here dot means differentiation with respect to t, v — x. 

Take the dimensionless variables rj and r, dimensionless force / and di- 
mensionless radiation parameter 7 and introduce scale multiplies a, b, p: 

2 

6 CL 

a ^ Xci, b ~ Xci, Xci = — ^, P = T 

52 2 

X = ar], ct = bT, f = F -, 7 = —Xd, 

amc^ 6b 

here x^i - the classical particle radius. 

In terms of new variables equation (1) becomes 



V 



(1 - (p?7)2)3/2 



/ + 7 



V 



+ 



1 — {pfjYY (1 — {pr]YY 



here dot means differentiation with respect to r. 
With the help of relativistic "velocity" u: 



u — 



V 



v/{cp) 



"acceleration" w: 



(1 - (p77)2)V2 (1 _ (^;/c)2)V2 

du 



w 



dr 



and dimensionless proper time s: 

c{l-{v/c)^y/^dt ^^^^^>> dr ds 
equation (2.2) can be put in more simple form: 

dw 



ds 



(2.2) 



(2.3) 
(2.4) 



(2.5) 
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The "solution" of (2.5) is obvious: 

w = w{s) = — exp(s/7) j dxf{x) exp{—x/'y) (2.6) 

^ so 

with "initial" value Wq: 

^0 = ( (!_ J^).)3/2 ) = ^ Jdxfi^) (-^/7) 

Integration of (2.4), taking into consideration (2.3) and (2.6), yields the 
following "solution" (strictly speaking, the integral equation) for particle 
velocity v: 

1 



In 



{l + v/c){l-vo/c) 



p ^ (l-^;/c)(l + Vc) 

s s 

j dzf{z) - exp (s/7) j dzf{z) exp {-z/'^) + 7^0 (exp (s/7) - 1) (2.7) 


here Vq- the "initial" velocity. 

The form of Lorentz-Dirac equation similar to (2.7) was given in [2,3,4,7]. 
Our form (2.7) is convenient for analysis - from it immediately follows that: 

(i) The peculiar features of Lorentz-Dirac equation do not qualitatively 
depend on scale multiplies a, 6, p, so these features are valid as for "small" , 
so for "large", "classical" distances. 

(ii) If conditions of a problem permit to consider the limit s ^ 00 (all 
integrals in (2.7) are not divergent for s < 00 ) then all "solutions" (2.7) 
must be "runaway" - \v\ — > c, with one exception. 

(iii) This exception is the particular case of zero asymptotic value of 
"acceleration" w: Sq = 00 in (2.6) and 

1 °° 

Wq = — dxf{x)exp(—x/^) (2.8) 
7 J 
' 

With (2.8), the R.H.S. of (2.7) takes the form 

s 00 00 

j dzf{z) + exp (s/7) j dzf{z) exp {-z/'j) - j dxf{x) exp {-x/-f) (2.9) 
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Then for s ^ oo and for "well defined" force / the velocity v does not 
reach c. But the price for this is the preacceleration and backward in time 
integration (see equation (2.9) ) with accompanying paradoxes (some of them 
are discussed in [7,8,9]) (see also section 3.). 

(iv) In literature the "solution" (2.9) is often called "physical", but from 
equation (2.7) it is easy to sec that (2.9) is unstable under small deviations 
of "acceleration" w from zero value at infinite "future": due to (2.7) these 
initially small at s = +oo deviations 5 grow at least as e^. 

The instability of the "physical" solution one can also verify with the help 
of numerical calculation: the instability was found in [10] and was checked 
by author for the force / formed by external charge at rest. 

Following (iv) one can state that there are no stable "nonrunaway" solu- 
tions of Lorentz-Dirac equation, at least in one-dimensional case. 

3. 

Here we present one more paradox of preacceleration: the formation of 
unavailable area of initial data, i.e. the formation of "event" horizon, absent 
in classical equation without radiation reaction. 

Consider for simplicity the nonrelativistic case (with units c—1). 

Let the point particle with mass m and charge e move (following the clas- 
sical nonrelativistic equation without radiation reaction) under the influence 
of the external force F along the x-axis: 

(fx _ F{x) 

Let's take the force F{x) in the form of a positive step (this form of the force 
is often used in the literature - for ex., [2,3]): 



0, X < —Xq ^ 

mA —Xq < X < 
0, 0<x I 



(3.2) 



here A > 0, > 0. In the method of backward integration one must take 
the "final" data (i.e. for t — > oo) and zero final acceleration. Thus for (3.2) 
we assume free motion of the particle in the future with velocity v: 

X — vt (3.3) 
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Suppose that the point x = is achieved at t = 0, and the point x = —Xq 
- at t = to < (the value of to we shall find below). 

Then the integration of (3.1-2,3) with appropriate boundary conditions 
(for the second order equation (3.1) the position and the velocity of particle 
must be continuous) yields: 



X 



f vt, t>0 \ 

vt + to<t<0 
\ ut + b, — oo < t < to J 



(3.4) 



Values of u, b, to are determined from the matching conditions: 

—v + u 



u 



j — ^ -l- ^ 

^yv'^-2xoA, to = — , b^-xo-uto (3.5) 



Following (3.4) in the regions free of force (a; > and x < Xq ) the particle 
motion is free. 

Consider the case of small velocity ^ — u at the point x — —xq: 



u 



= yv2^^2^=e, e^O (3.6) 



We see that though the value of u is small, the particle can reach from the 
future all points on the x-axis: —oo < x < +oo. 

Consider now the same particle with the same final data moving under 
the same external force F (3.2) but obeying the nonrelativistic Lorentz-Dirac 
equation 

(fx _ (fx _ F{x) 
dt^ "dt^- m 

here k {k ^ x^-the classical radius of a particle. A; > 0) and the second term 
on L.H.S. of (3.7) deal with the radiation force. 

As the equation (3.7) is of the third order, the acceleration of the particle 
also must be continuous. Then the solution of the above problem for equa- 
tion (3.7) with appropriate boundary conditions yields (the point x = is 
achieved at t = 0, and the point x — —Xq - at t = ti < ) 



x^\{v + kA)t + At^/2 + k^A{l - exp {t/k)) 
ut + b + cexp {t/k), 



t > 
ti < t < 
t < ti 



(3.8) 
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Values of u, b, c, ti are determined from the matching conditions: 



—Xo — uti + b + cexp{ti/k) (3.9a) 

-xo = (f + kA)ti + A{tif/2 + PA{1 - exp (ti/k)) (3.96) 

v + kA + Ati~ kA cxp (ti/k)) = u + (c/k) cxp (h/k) (3.9c) 

A — A exp (ti//c) = (c/A;^) exp (h/k) (3.9d) 



Term with exponent in (3.8) for x < —xq (i.e. "before" the action of the 
force) describes the effect of preacceleration. 

From equation (3.9d) immediately follows that c > for ti < 0. 

Consider the case of zero particle velocity ^ at the point t — ti <0, x — 
-Xq: 

(IT 

-^^u + {c/k) exp {tr/k) = (3.10) 

In (3.10) the vahic of u must be negative because the value of c is positive. 
Equation (3.10) with the help of the system (3.9) can be rewritten as 

^V2-1-P = exp(^)(^-1) (3.11) 

here z — ti/k < 0, p — xo/{k'^A) > 0. For our goal it is important to note 
that the equation (3.11) always has solution for z < and this solution varies 
weakly with small changes in parameters of the problem under consideration. 

Consequently if we consider the similar to (3.6) case of small particle 
velocity at the point x — —xq: 

^{h)^e, e^O (3.12) 
then t2 ~ ti and, as for (3.10), 

u <0, c> 

. These inequalities lead to the conclusion that the particle velocity in the 
free of force region {x < —Xq) - ^ = u + {c/k) exp {t/k), with positive value 
€ at X — —xo, must inevitably take (with time decrease) zero value at some 
moment t = tmin < h and some point x — Xmin < —xo- 

So all X on the left of Xmin' x < Xmin, become unavailable on contrary 
to solution (3.6) of the equation (3.1), where all x-axis is available for the 
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particle. Thus the event horizon is formed: not all initial data, physically 
possible from the point of view (3.1), can be achieved by backward in time 
integration method. 

This paradox (as all the other) makes unrealizable the desire to solve 
problems of Lorentz-Dirac equation using preacceleration and backward in 
time integration. 



If to accept that there are problems with the Lorentz-Dirac equation, 
then one can try to solve them considering radiation forces of another type. 

In the literature there are many examples of different forms of radiation 
forces. Let's mention some of them. 

4.1 

Following the works of Tcitelboim school [5,6], when a charge is acceler- 
ated by a external field, it has a bound momentum given by 

fly} 

= P;„,, + Pl,^, = mu'-m-k-—, 2eV(3m) 

With it the Lorentz-Dirac eq. must be rewritten as 

dP' ■ ,du'' dur,. 

F"- + m- ku" 



ds ds ds 



4.2 

In works of Plass [3] and Rohrlich [1] was shown that taking into ac- 
count the asymptotic conditions (zero values for acceleration at infinity) the 
Lorentz-Dirac eq. can be rewritten in the form of integrodifferential (nonlo- 
cal) equation: 

dv} f 
m——{s) = I daK\s + k ■ a)e~°' 
ds J 


d^v"' 

K'^F'-m- ku'un-r^ 

ds^ 
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4.3 



One can rewrite the Lorentz-Dirac eq. using the iteration procedure con- 
sidering the radiation term as small one (for recent paper see, for ex., [20]): 

m— = F'; 
as 

m—— ^ F' + m- k(— u un—;—) = -Fi 

as as as 

m— = F' + m- ki—^ - u'un—f^) = 
as as as 

...etc 



4.4 

One can try to consider general form of the eq. of motion for a point 
particle using the distribution theory (Lozada, [14]) 

with T*^- the total energy-momentum tensor for a point particle in an external 

electromagnetic field. 

From this general eq. under some assumptions one can derive in particular 
the Bonnor's theory [21] with total nonconstant mass of a particle 

m — m(s); m—— — 
as 

dm , du" dur, 

—— — m ■ k- 



ds ds ds 

4.5 

One can draw the quantum field theory to these problems and to consider 
a quantum particle coupled to a quantum-mechanical heat bath (more, one 
can find the opinion in the literature that only quantum field approach can 
solve the problems of radiation reaction. I do no share this point of view 
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as the only one possible). For such a system the description in terms of 
the quantum Langevin equation has a broad application (a far as I know, 
brownian model of synchrotronically radiating particle was first suggested 
by Sokolov and Ternov [18], for recent developments see , for ex.. Ford et 
al [18]). Then in the nonrelativistic case the macroscopic eq., describing the 
particle motion, is 

t 

mx+ J dt'n{t-t')x{t')^ f + Fe^t 

— oo 

here jj, - the memory function, / - a random force with mean zero. 
For microscopic eq., see, for ex., Efremov [18]. 

All the above approaches use the retarded form of self-interaction. What 
new can give the consideration of a superposition of retarded and advanced 
self-interactions? For this sake let's take the model of extended (in some 
sense) particle. 



4.6 



Consider the hydrodynamic model of an extended particle [15]. 
Then the particle is described by the mass density m ■ f{t,x), charge 
density p{t,x) and current density j(t,x), obeying the continuity equations 

df ^ d{v ■ f) 
at ox 

here v = v{t,x) is the hydrodynamic velocity of moving extended particle. 

Let the particle move under the external force F(t, x) along x-axis. Then 
the relativistic equation of its motion reads (we choose the units c—1): 

m J dxf{t,x) ^^ + v^^u = Jdxp{t,x)E{t,x)+Jdxf{t,x)F{t,x) (4.2) 

here v — v{t,x), u — u{t,x) = v/\/l — v"^, E{t,x) -is the electric field, 
produced by moving particle (the Lorentz force is absent in one-dimension 
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case under consideration and internal forces give zero contribution to total 
force) : 

and electromagnetic potentials and A are 



J \x — x'\ 



A{t,x) = j dx'dt'j^^AaSi + 652), (4.4) 
with retarded and advanced delta-functions 

5i = 5{t' -t+\x-x'\), ^2 = 5{t' -t-\x-x'\) 

and a,b - constants. 

(The role of retarded and advanced interactions in electrodynamics is 
vividly described in the textbook [16].) 

Substitution of (4.4) in (4.3) and integration by parts with the help of eq. 
(4.1) (taking zero values for integrals of exact integrands in x', i.e. / dx'^{p- 
••) = 0), yields 

E{t,x) = I -^^^^(p{t\x')^^^{a6i + b62)+j{t\x'){a6i-b62)] 

J \X X \ \X X J 

(4.5) 

Similar integration by parts for LHS of (4.2) gives the common result 

dP r 
LHS = — , P = P{t) ^mj dxf{t, x)u{t, x) (4.6) 

here P - the particle momentum. Thus the eq. of motion reads 

Pself ~l~ Fexti 

F.eif^ J dxp{t,x)E{t,x), Fe,t = J dxf{t,x)F{t,x) (4.7) 

This eq. of motion has no second derivative of particle velocity; also there is 
no need in mass renormalization. 
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If the extended particle is compact, we can use in (4.5,7) the standard 
expansion in powers of \x — x'\ (see, for ex., [2]): 

with e = ±1. 

Thus in nonrelativistic case we get the known resuh: 

^ , /■ dxdx' , ,, , ,dv(t.x') 

Fseif = -{a + b) J ^^_^^^ p{t,x')p{t,x) ^ + 

^(a - 6) y dxdx'p{t,x')p{t,x)^^^i (4.8) 

The first term in (4.8) is considered in hterature as "-(electro dynamic field 
mass) X (acceleration)", and the second - as radiation reaction force. 
The total change in particle momentum is 

r dP 

AP = P(oo)-P(-oo) = y 



dt- , 

dt 



and the change in particle momentum due to self-interaction is 

^Pself = / dtFself 

Thus 

AP = APseif + I dtF,^t (4.9) 
Substitution of (4.5) into (4.7,9) gives 



^ = F^elf= fdt'dxdx'/^. 

dt ^ J x-x'^ 



\p{t',x')-—-{aS, + hd^) +j{t',x'){aSi - hd^) (4.10) 

APseif = j dt\RHS of (4.10)] (4.11) 
The solution of (4.1) we can write as 
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Then integration by parts in (4.10,11) with the help of (4.12) gives the fol- 
lowing result: 



AP,eif = J dtdt'dxdx'^{t, x)^{t', x' 



X — x' 
\x — x'l^ 



-6(a5, + b6,) + 6\x-x I ^a— + b—j-2\x-x \ (^a^ + 

(4.13) 

In (4.13) the integrand is antisymmetric under transformations 
if 

a = 6 (= 1/2) 

. Then the whole integral (4.13) has identically zero value: 

^Pself = 

So for an extended in one dimension particle the total change in particle 
momentum due to its self-interaction is zero (if is taken the half-sum of 

retarded and advanced interactions). 

For a = b the radiation term in Fg^if (4.8) is identically zero. 

Similar result holds for energy balance: from equations (4.1,2) one can 
derive 

-Wkin^ A,elf + Ae^t (4.14) 



dt 



with 



m ■ f 



Aseif = j dx{j ■ E), Aext = j dx{Fe^t ■ v) 



and 



^Wkin = I dtjWkin = J dtAself + j dtAkin (4.15) 

Substitution of (4.12) into (4.14) gives for / dtAgeif identically zero result in 
the case 

a = 6 = 1/2 
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Then the natural question arises - where is the source of radiated energy 
in this case? - The answer is obvious and comes from eq. (4.9, 15) - the 
source of radiated electromagnetic energy is the external force work. 

If a 7^ 6 then APgeif is not zero and more, the sign of APgeif is not 
identically negative - for some processes it can be negative and for the other 
- not. The latter case one can consider as " antidamping" . (For example, 
if p and j have the form of moving "extended rigid" body with velocity 
v{t) = dR{t)/dt: p{t,x) = Aexp[-a{x - R{t)y], j = v{t)p{t,x), A,a- 
const, then for a = 1, 6 = (retarded self-interaction) APgeif > if v > 
for all t.) 

Consequently we see that the problems of radiation reaction are connected 
not with the chosen form of radiation force but rather with the form of chosen 
self- interaction - retarded or/and advanced. 

4.7 

Let's postulate in section 4.6 the following relations 

f{t,x) = j dvf{t,x,v) 

p{t,x) = Q J dvf{t,x,v) 

jit, x) = Q J dw fit, X, v) (4.16) 

with f{t,x,v)- the distribution function of extended in " f -dimension" par- 
ticle. Thus we introduce the stochastic character of interaction of particle 
with its self-electromagnetic field ( with "heat bath" in spirit of works [18, 
19]). Distribution function f{t,x,v) must obey the continuity equation: 

df{t,x,v) ^ d{vf{t,x,v)) ^ d{vf{t,x,v)) ^ ^ 
dt dx dv 

With it the equations of sec. 4. 6, describing the system "particle + self-field", 
take the self-consistent nonlinear form [19] and the results of section 4.6 
remain valid with the following interpretation: 

the stochastic self -interaction of particle can give zero contribution to total 
changes in particle energy and momentum. 

5 



14 



In the literature one can often find the statement that the insert of Egeif 
(or Fseif, Wem,seif = / dV {E^ + H"^) sei f / , etc.) into equation of particle mo- 
tion (4.2,7) is an obvious procedure justifying by the principle of extremum 
of action. Meanwhile the latter is formulated strictly only for closed systems 
or for systems, interacting with known external sources possessing known 
equations of motion. Systems with dissipation (with "damping" as a con- 
sequence of emission to infinity some amount of energy) may not obey the 
principle of extremum of action. Consequently the choice of the forms of 
Eseif (or Fseif, Wem,seif, etc), dealing with dissipation, can be considered as 
some supplementory hypothesis. 

Consider the equations (4.14) and (4.9) together with the balance equa- 
tion of classical electrodynamics: 

^M^em = -/-^d/ (4.14) 

where We„, = / dV{E'^ + i7^)/87r - the total energy (with the self-energy) in 
the volume V with surface S containing the moving source of electromagnetic 
field, and / - is the fiux of Poynting-Umov vector through the surface S. 
Then the latter can be put in the form 

^Wtot = -/ + Ae^t 
at 

with Wtot = Wem + Wkin Or 

^Wtot ^-Jdtl + J dtAe^ (5.1) 

This equation is verified in experiments. More precisely, the R.H.S.'s of this 
equations - because the explicit value for Wtot is unknown in the case of point 
charged particles. Nor is known in the experiments the explicit equation of 
motion of radiating charged particle. So one can take for self-interaction the 
superposition of retarded and advanced interactions without violating the 
main equation (5.1). 

Thus we can reformulate the hypothesis of section 4 in the following form: 
For an extended (in some sense) particle the contribution of particle self- 
interaction to the total energy-momentum balance of system " particle -h field" 
depends not on the choice of the form of radiation force, but rather on 
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the choice of superposition of retarded and advanced fields, describing self- 
interaction of an particle. 

Then in classical electrodynamics (in the classical theories of scalar, grav- 
itational, etc. fields) the use in eq. of particle motion radiation reaction force 
in Lorentz-Dirac or another form is not the inevitable procedure, dealing with 
the fact of radiation, the radiation of field energy can be possible due to the 
work of some external force with zero value for '' radiation force" . 

In favor of this hypothesis speaks the following example from the classi- 
cal theory of a scalar field. Consider the model of a scalar field with self- 
interaction in the form: 

L = —dp(pd^(j) + g ■ n ■ (f) — m ■ n, dp{nu^) = 

OTT 

Consider the self-expansion of a scalar-charged sphere. Then the eq. of its 
motion is determined by the discontinuities of the total energy-momentum 
tensor across the sphere. This tensor is constructed in rather simple way 
due to the spherical symmetry of the problem. Then the eq. of motion 
is self- consistent - there is no need to add to it radiation force - radiation 
effect is automatically taken into account if one knows the total tensor on 
each side of the sphere (it is necessary to note that, contrary to classical 
electrodynamics, the expanding scalar-charged sphere in the classical scalar 
field theory do radiate). The numerical analysis of such expansion shows [22] 
that there are the solutions that one can interpret as "runaway": t; — > ±1. 

Thus the existence of such solutions does not bound up with the forms 
of radiation force. 
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